Basins of attraction of a nonlinear nanomechanical resonator 
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We present an experiment that systematically probes the basins of attraction of two fixed points 
of a nonlinear nanomechanical resonator and maps them out with high resolution. We observe a 
separatrix which progressively alters shape for varying drive strength and changes the relative areas 
of the two basins of attraction. The observed separatrix is blurred due to ambient fluctuations, 
including residual noise in the drive system, which cause uncertainty in the preparation of an initial 
state close to the separatrix. We find a good agreement between the experimentally mapped and 
theoretically calculated basins of attraction. 



In the last few years the dimensions of mechanical 
devices have been scaled deep into the submicrometer 
regime. This decrease in size has resulted not only in the 
increased detection sensitivity of extremely small physi- 
cal quantities, such as zeptogram-scale mass [l[ and single 
electron spin , but also in an enhancement of the sig- 
nificance of nonlinear dynamics in such devices 0]. The 
growing significance of nonlinearities in high-frequency 
nanomechanical devices has consequences not only for 
fundamental studies of nonlinear dynamics, but also for 
advances in sensing applications. It was recently shown 
that precision of some experimental measurements on 
nanoscale can be improved by deliberately operating the 
system in the nonlinear regime. For example, a non- 
linear resonator can be employed to suppress amplifier 
noise in an oscillator circuit [4], noise- induced switching 
between two stable states in a nonlinear beam resonator 
enables precision measurement of the resonant frequency 
, and the sensitivity of a resonator for mass detection 
can be improved when the resonator is driven into a re- 
gion of nonlinear oscillations ^ . Finally, in a Josephson 
junction, which is dynamically similar to a mechanical 
resonator in nonlinear regime, the bistable state of the 
nonlinear system can be used as a bifurcation amplifier 
to perform a non-dissipative, low-back-action measure- 
ment of the phase across the junction [7]. When nanome- 
chanical devices reach the quantum- limited regime 8], 
a nanomechanical version of such an amplifier could be 
used for a similar sensitive low-back-action measurement 
of the state of a quantum mechanical resonator. Non- 
linear response of nanomechanical resonators could also 
be used to detect transition from classical to quantum 
regime 0, It is therefore important to understand 
nonlinear dynamics of these systems well, so that we can 
fully realize their potential in expanding our experimen- 
tal capabilities. 



Although some work has been done with parametric 
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^ [, the majority of nonlinear nanoscale sys- 
tems that have been studied are directly driven 0, 13, 14 1 
and the dominant nonlinearity in the restoring force is cu- 
bic, also known as Duffing nonlinearity. When a system 
is driven strongly, the Dufhng nonlinearity causes the 
resonance response curve to become asymmetric. The 
resonance is pulled either to the right for positive, also 
known as hardening, nonlinearity (e.g. geometric non- 
linearity [l^. Figure [Ijb)) or to the left for negative, 
or softening, nonlinearity ( e.g. nonlinearities of material 
[iHl, capacitive, or inertial [l3| origins). When the reso- 
nance is pulled far enough to one side, hysteretic behav- 
ior is observed as two stable states appear in the system 
[ist . The stable states, known as "attractors" or "fixed 
points" , correspond to the points in state space to which 
the system converges with time. For each attractor, a set 
of initial states that dynamically evolves to that attrac- 
tor forms its basin of attraction, which is separated from 
the rest of the state space by the separatrix curve. The 
dependence of fixed points and basins of attraction on 
driving frequency and amplitude has been exploited for 
precision measurement applications mentioned earlier. 

There have been very few experimental studies of 
basins of attraction because following the evolution of 
initial conditions in low-frequency macroscopic systems 
is usually very time-consuming and system parameters 
tend to drift over the course of many data-taking runs. 
Previous mappings of basins of attraction [l^j EO] used 
the method of stochastic interrogation, where the system 
is stochastically perturbed and initial states are sampled 
at random without fully covering the basins. 

We demonstrate the ability to systematically prepare a 
nonlinear, Duffing-type nanomechanical resonator in the 
required set of initial states and map the basins of at- 
traction of its two fixed points with high resolution. Our 
experiment is enabled by the fact that the relevant time 
scale per data point, ~ Q/ fo, is very short for very-high- 
frequency nanometer-scale devices. This allows us to take 
many data points in a relatively short time with mini- 
mal drift in the parameters of the system. We also ob- 
serve that the separatrix changes shape for varying drive 
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FIG. 1: (a) Experimental layout. The initial drive prepares 
an initial state of the platinum nanowire resonator (shown in 
the SEM photo), a 5 ns RF switch is then flipped to connect 
to the final drive, and the state of the device is measured by 
a lock-in amplifier after mixing down to a low (1 kHz) fre- 
quency and filtering the residual RF signal. 

(b) Vibration amplitude versus frequency, for various driving 
powers at the sample (-90, -85, -80, -75, -70 dBm, or nor- 
malized to the critical drive 14=28.4 ^V, l^/K=0.249, 0.443, 
0.788, 1.401, 2.492) showing the onset of nonlinearity in the 
platinum nanowire resonator. We plot the response normal- 
ized to the critical amplitude, ac, versus normalized detuning 
frequency a = 2Q{f/fo — 1). The backbone curve (dashed 
line) connects the maxima of the resonance curves and fol- 
lows {ap/acY ~ \/icT/2. 

(c) The state of the nanowire resonator is first prepared in an 
initial state on the initial response curve (grey) by choosing an 
appropriate drive strength for the fixed detuning frequency, 
Od = 4.26. After the RF switch is fiipped to connect the final 
drive, the response curve changes to the one shown in black 
and the initial state evolves to either the high-amplitude state 
or low-amplitude state. The small bias voltage, applied to the 
nanowire together with the final drive, is chosen so that the 
hysteresis loop of the initial curve is at higher frequency than 
the operating frequency. 



strength, so that one of the basins becomes progressively 
smaller and eventually disappears. The mapped basins 
of attraction show good agreement with theory. How- 
ever, the observed separatrix is blurred due to ambient 
fluctuations, including residual noise in the drive system, 
which cause uncertainty in the preparation of an initial 
state. 

The device used for mapping the basins of attraction, 
a doubly-clamped platinum nanowire, is shown in the 
scanning electron microscope (SEM) photograph in Fig- 
ure [T]J a). The nanowire, with a length L of 2.25 |J.m and 
a diameter of 35 nm, is grown by electrodeposition of 
platinum into a nanoporous membrane [211] . Gold con- 
tact pads on both ends and a gate are fabricated using 
electron beam lithography, and about 150 nm of the sub- 
strate is subsequently etched away in hydrofluoric acid to 



suspend the device [l3| . We actuate and detect the vibra- 
tion of the nanowire magnetomotively [2^ in a magnetic 
field, i? = 8 T, in a cooled (to about 20 K) probe in vac- 
uum . The magnetic field is applied perpendicular to the 
device so that the vibration is in the plane of the gate 
electrode. At low driving powers the resonance curve is 
linear, and we extract a resonant frequency of 45.35 MHz 
and mechanical quality factor of 6045. The resonant fre- 
quency is higher than the expected 17.71 MHz for this de- 
vice geometry, most likely due to differential thermal con- 
traction between the silicon wafer and the gold contacts 
that results in residual tension. The ratio of electrome- 
chanical impedance, Rem, to electrical impedance, i?e, 
is 0.222, which indicates the presence of significant eddy 
current damping [2^ . The quality factor corrected for the 
eddy current damping is Qq = (5/(1 — Rem I Re) = 7770. 

The resonant response of the nanowire to different 
drives is shown in Figure \V[h) . With increasing drive 
power, the resonance is pulled to higher frequencies at 
large amplitudes, ultimately forming a hysteretic region. 
The nonlinearity of the device is fully characterized by 
the critical amplitude the point where the resonance 
curve develops infinite slope, |a=a<,= —oo. The the- 
oretical curves (thin solid lines) in Figure mb) are gen- 
erated using the critical amplitude, 

(1(2 ; the only fit- 
ting parameter. We determine Oc by fitting the back- 
bone curve that connects the peaks of resonant curves for 
different drives to the theoretical expression (op/cc)^ = 
\/3cr/2 [ill, where a = 2Q(///o - 1) is the detuning 
frequency scaled by the width of the resonance and Up 
is the peak amplitude. The critical amplitude value 
flc — 2.68 nm, extracted from the experimental data in 
this manner, is in reasonable agreement with the value 
calculated for our nanowire geometry Q when the round- 
trip loss in the experimental circuit is taken into account. 

The dominant source of nonlinearity in doubly- 
clamped NEMS resonators is the additional tension in 
the beam that appears when vibrations are sufficiently 
large. This extra tension gives rise to a cubic nonlinear- 
in the spring constant term in the equation of motion 
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iving it a Duffing oscillator form: 

X + '^x + ljq{x + ax^) = Fcos{ilt). 



(1) 



Here, x{t) is the displacement of the beam, uq ~ 27r/o is 
the resonance frequency, a = 2\/3/(9a^(5) is the nonlin- 
earity parameter [2J], 17 = 27r/ = {a /{2Q) + 1)loq is the 
driving frequency, and F is the force per unit mass act- 
ing on the resonator of mass m. The driving force in the 
magnetomotive transduction scheme is the Lorentz force 
that acts on the nanowire when a current Id{t) is passed 
through it in a magnetic field, F = LBId{t)/m. When 
the driving voltage amplitude is Vd, the driving current 
is given by Id^Vd/Re- 

As illustrated in Figure [Ija) , we prepare the system 
by exciting it with an initial drive amplitude, Vi, and 
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then switching to a final drive ampHtude, Vf. Two RF 
sources (HP 8648B for initial and SR DS345 with a fre- 
quency doubler for final drive) are tuned to the same 
fixed frequency off resonance, — 4.26. Their internal 
clocks and that of the local oscillator (LO) are synchro- 
nized with their 10 MHz clock reference. The phase of 
the final drive lags behind the initial drive phase by the 
phase difference 0. By changing (/) while holding the ini- 
tial drive amplitude, Vi, constant, we can prepare the 
resonator in the initial states corresponding to a circle 
in state space. By also stepping the initial drive values, 
Vi, we can cover a disk of initial states in state space. 
We switch rapidly from the initial to final drive using a 5 
ns (< l//o ^22 ns) RF switch (Mini-Circuits ZASWA-2- 
50DR). After the switching occurs, we measure the final 
amplitude of the oscillator and mark it as a low or a high 
final amplitude. In order to access a continuum of ini- 
tial states, we apply a small DC bias voltage of Vdc ~ 10 
mV to the wire in the final state. The capacitive inter- 
action with the gate lowers the resonant frequency of the 
final state so that the hysteretic frequency response 
of the final state occurs at the same frequency as the 
single- valued resonant response of the initial state (Fig- 
ure [ijc)). Without this technique, the resonator cannot 
be prepared in the initial states corresponding to the un- 
stable branch of the initial drive resonance and an annu- 
lus of these states in state space would not be accessible 
in the experiment. 

To map the basins of attraction in this manner, the 
initial states were driven with —90 to —50 dBm in 60 
concentric circles with 60 phase points per circle, corre- 
sponding to a displacement range of to 2.5ac, where 
ttc is the critical amplitude extracted from experimental 
data as shown in Figure [T](b) . Each of the initial states 
was marked according to the attractor that it evolved to 
after the switch was flipped: blue for the high-amplitude 
state and yellow for the low-amplitude state. This data 
was re-rastered using a nearest neighbor search algorithm 
to create the continuous color plots shown in Figure [21 

For very low final drives, there is only one state the 
resonator can occupy (refer to Figure [TJb)). As the final 
drive starts exciting a nonlinear response, a second basin 
representing the high-amplitude stable state appears, but 
subtends a smaller fraction of the state space (Figure 
[2^1)). When the final drive amplitude is increased, the 
high-amplitude basin grows, and the low-amplitude basin 
gets progressively smaller and eventually disappears (Fig- 
ure [2l(6)). The disappearance of the low-amplitude state 
can also be observed in Figure [Ijb): a large final drive re- 
sults in a wider hysteretic region that moves to the right, 
where the low amplitude branch starts at cr > cr^, so 
only the high-amplitude state is available to the system 
operating at a^. 

In order to calculate the location of the fixed points 
and the separatrix in state space, we obtain the time- 
dependent solution x{t) that describes the response of 
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FIG. 2: Basins of attraction of a nanowire resonator at a 
fixed frequency, Od ~ 4.26, for increasing final drive values, 
V>/V'c=(l)1.867, (2)2.049, (3)2.237, (4)2.434, (5)2.640, (6)2.741. 
Blue and yellow colors indicate the final high- and low- 
amplitude states, respectively. The data consists of 60 
concentric circles with 60 points each, corresponding to a 
displacement range of to 2.5ac. This data is converted into 
a continuous plot using a nearest neighbor search algorithm 
to fill out the space between data points. Theoretical fixed 
points, saddle point, and separatrix curve are indicated 
by black points, black cross, and the dashed black curve 
respectively. 



the nonlinear system to a disturbance. In our devices, 
we can separate the dynamics described by the equation 
of motion ^ into two parts: the fast dynamics on a time 
scale of l/wg, corresponding to the fast oscillations of the 
undamped harmonic version of the system; and the slow 
dynamics on a much longer time scale Q/wq, associated 
with a slight detuning a of the driving frequency from 
resonance as well as damping and nonlinearity (method 
of multiple scales [3]). Then the solution to the equation 
of motion ^ can be written as xo{t,T) = ^(T)e''^°* -I- 
yl(T)e~*""*, where time variable t characterizes the fast 
dynamics and T, the slow dynamics. The slowly varying 
amplitude A{T) = {X{T) -)- ir(T)) exp(iwocr/2T) obeys 
the envelope equations: 
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We have assumed here that l/Q ^ 1, and the slowly- 
varying amplitude approximation implies that A{t) terms 
are negligible compared to A terms. 

The slowly-varying amplitude equations describe the 
nonlinear dynamics of the system and allow us to deter- 
mine the location of two attracting fixed points and one 
metastable saddle point in state space for different values 
of parameters F and ct, which can be extracted from the 
experimental data. A set of points in state space that 
evolve into the saddle point defines the separatrix. To 
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calculate the separatrix, we evolve the initial conditions, 
lying close to this fixed point and along the negative- 
eigenvalue eigenvector (which is obtained by linearizing 
the above equations around the saddle fixed point), back- 
wards in time according to the full equations The 
curves generated by this procedure constitute the sepa- 
ratrix. To compare these theoretical calculations to the 
experimental data in Figure [21 we scale computed ampli- 
tudes by the value of the critical amplitude, Oc, calculated 
from equations 

The black points in Figure [2] are the theoretical fixed 
points, the black cross is the saddle point, and the dashed 
black curve corresponds to the theoretical calculation 
of the separatrix for the experimental parameters used: 
ad — 4.26, {a/ac)max = 2.5, and the first five final drive 
values scaled by critical drive, Vf /Vc, listed in the cap- 
tion. We find good agreement between the experimental 
data and our theoretical calculations for the basins of 
attraction. 

The separatrix observed in the experiment is blurred 
due to environmental noise affecting the system. The 
following analysis implies that most of this noise results 
from the voltage noise in the drive circuit. These fluc- 
tuations perturb the initial state of the resonator and so 
cause uncertainty in the preparation of the initial state. 
For the states near the separatrix, switching the system 
to the bistable regime and forcing the initial state to 
project on either of the two stable states amplifies this 
jitter and results in the noisy separatrix in Figure [21 In 
principle, the full solution for the probability distribu- 
tion of the initial state of the Duffing resonator with noise 
would be obtained by solving the Fokker-Plank equation. 
Since the system is operated away from the bifurcation, 
the time scale on which noise affects the system is much 
slower than the ringdown time of the resonator and the 
dynamics can be analyzed in the small-noise approxima- 
tion. We can then consider two scenarios for how exter- 
nal noise can affect the system, either additively or as 
parameter noise. 

If noise in equations ^ for the initial state is additive, 
the evolution of the initial condition follows the deter- 
ministic path perturbed by small jitter until the system 
reaches the neighborhood of an attracting fixed point, 
corresponding to the initial state. The action of this noise 
is to perturb the system around the fixed point within a 
characteristic radius 5a. If the additive noise in equation 
(jl]) is white with power spectral density 7, then the effec- 
tive temperature of the noise in the system described by 
equations ([2]) is g^^j ■ From the equipartition theorem, 
5a is given by 
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We can estimate the uncertainty in the preparation of 
the initial amplitude from the observed blurring of the 



separatrix, 5a/ ac, in Figure[2]to be about 10%. The force 
noise that would result in such blurring is 5 pN/VHz. For 
the Lorentz force due to the magnetomotive drive, this 
force noise translates to a voltage noise of 18 nV/-\/Hz. 

Alternatively, the noise can be present in the param- 
eters of the system. The noise in the magnitude of the 
drive, F, results in the fiuctuations of the position of the 
initial fixed point. Near the origin, where a/oc is small, 
the variation in the drive voltage 5V/Vc that would cause 
this variation in amplitude is calculated from equations 
(HI) to be 5V/Vc = (V3/2)a/ac. For the estimated blur- 
ring of the separatrix, 5a/ ac, of about 10%, the variation 
in the drive amplitude is 5V/Vc = 8.7%. (Farther away 
from the origin, up to a/ac ~ 2.5, fiuctuations in the 
drive voltage have less effect on ffuctuations in the am- 
plitude due to nonlinear suppression.) Therefore, the ob- 
served blurring translates to fiuctuations of i5F « 2.5 /xV 
for the measured critical drive, Vc = 28 /iV. The relevant 
noise bandwidth for this resonator is irfo/Q — 24 kHz. 
The noise spectrum that would account for the 10% fluc- 
tuation in amplitude is then 16 nV/\/Hz. This result 
is consistent with the above estimate of fiuctuations due 
to additive noise. The residual voltage noise from the 
initial-drive function generator and the rest of the drive 
circuit of about 5 /iV/VHz is attenuated by 51 dBm by 
the RF switch and additional attenuators (not shown), 
and results in a voltage noise of 14 nV/VHz at the sam- 
ple. The drive-circuit noise therefore accounts for most 
of the observed fluctuations near the separatrix. 

The presence of noise can also cause transitions from 
one fixed point to another 2^, 27 1 . We indeed observe the 
same noise-induced switching between two stable states 
as in References Q and [2^, where the noise from the 
drive source has an effect of shrinking the size of the 
hysteresis loop and inducing transitions of a resonator 
from one state to the other near the bifurcation points. 
By adding noise to the resonator drive and recording the 
statistics of the time it takes for the system to switch 
when it is near the bifurcation point, we find that the 
transition rate varies as exjp{—Ea/i'), where v = ksTeff 
is the noise power and Ea is the height of the energy bar- 
rier that the system needs to overcome for the transition 
to happen. The energy barrier depends on the distance to 
the bifurcation point 14: Ea ^ [V — VbY ■ We have mea- 
sured the critical exponent 5 to be 1.8±0.3, which is close 
to the theoretically predicted value of 3/2 in the region 
we operate in [2^, |23, [2^. Measurement of transitions 
induced by noise in the bistable regime of a nanoscale 
resonator could thus enable a very sensitive experimen- 
tal technique to detect mechanical fluctuations. 

The experimental mapping of basins of attraction of 
a nanowire mechanical resonator presented here fills a 
large gap in our understanding of nonlinear dynamics of 
nanoscale systems. Since the nonlinear regime is readily 
accessible in nanoscale devices, the details of dynamical 
behavior are now increasingly important for proper en- 
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gineering and analysis of these systems. The knowledge 
of basin dynamics and noise-induced transitions should 
also prove useful for precision measurement applications, 
such as nanomechanical bifurcation amplifiers 1121 , de- 
tection of transition to the quantum regime llC| . or a 
sensitive monitoring of intrinsic device noise processes. 
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